N=64;
j=sqrt(-1);
a=ones(N,1);
row=0;
for i1=1:2:N;
	a(i1) = sqrt(row+i1-1+.77777);	
	a(i1+1) = (row+i1-1+.22222)*(row+i1-1+.22222) / N - N/2;	
end;
a(1)=N+3;
a(3)=1-N;
A2=fft(a);
A=A2(1:N/2);
A(1)=A(1)+j*A2(N/2+1);
a1=reshape(a,N/2,2);
if N<=32
	A
	a1
end
max(abs(A-atrans))./(max(abs(A))/N)
max(abs(a-aitrans))./(max(abs(a))/N)
